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ABSTRACT 

An important class of formation theories for hot Jupiters involves the excitation of extreme orbital eccentricity 
(e = 0.99 or even larger) followed by tidal dissipation at periastron passage that eventually circularizes the plan- 
etary orbit at a period less than 10 days. In a steady state, this mechanism requires the existence of a significant 
population of super-eccentric (e > 0.9) migrating Jupiters with long orbital periods and periastron distances of 
only a few stellar radii. For these super-eccentric planets, the periastron is fixed due to conservation of orbital 
angular momentum and the energy dissipated per orbit is constant, implying that the rate of change in semi- 
major axis fl is fl' C!C a'/- and consequently the number distribution satisfies dJV/d\oga cx a'/-. If this formation 
process produces most hot Jupiters, Kepler should detect several super-eccentric migrating progenitors of hot 
Jupiters, allowing for a test of high-eccentricity migration scenarios. 
Subject headings: extra-solar planets - tidal friction 



1. INTRODUCTION 

The origin of gas-giant planets with orbital periods of only 
a few days - the hot Jupiters - is not understood. One hypoth- 
esis involves the following sequence of events: (i) the planets 
form at a few AU from their host stars, in approximately circu- 
lar orbits; (ii) some mechanism excites their orbits to extreme 
eccentricities (1 -e < 0.01); (iii) tidal dissipation during suc- 
cessive periastron passages removes enough orbital energy so 
that the planet migrates a factor of ^ 100 in semi-major axis, 
finally settling into a circular orbit close to the host star. Pos- 
sible excitation mechanisms include Kozai-Lidov (KL) oscil- 
lations (Wu & Murray 2003; Fabrycky & Tremaine 2007), 
planet-planet scattering (Rasio & Ford 1996; Nagasawa et al. 
2008), resonant capture during migration in multi-planet sys- 
tems (Yu & Tremaine 2001), and weak resonant orbital inter- 
actions (called "secular chaos" by Wu & Lithwick 201 1). We 
shall refer to these as high-eccentricity migration (HEM) sce- 
narios; they are of particular interest because they naturally 
predict frequent misalignment of the stellar spin and plane- 
tary orbit, consistent with recent observations of the Rossiter- 
McLaughHn effect (Winn et al. 2010). 

Most hot Jupiters have relatively small eccentricities: 60% 
of the known planets with orbital period P < 1 d and M sin / > 
0.25M/ (Jupiter masses) have eccentricities consistent with 
zero, and 90% have eccentricity e < 0. 1 . If the hot Jupiters are 
formed through the HEM process, this result implies that the 
timescale for the decay of the eccentricity from near unity to 
near zero is short compared to the age of the Galaxy. Since the 
star-formation rate is approximately constant over the age of 
the Galaxy, the distribution of migrating Jupiters with moder- 
ate or high eccentricity should therefore be in an approximate 
steady state. Moreover, for these eccentricities the energy dis- 
sipation per orbit is independent of eccentricity, since the dis- 
sipation occurs only near periastron and in this region all mod- 
erate and high-eccentricity orbits look like parabolae. There- 
fore, we can predict the eccentricity and semi-major axis dis- 
tribution in HEM independent of the details of the dissipation 
process. We do this in §2 and find that HEM requires the 
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presence of a large population of Jupiters with eccentricity 
e > 0.9, which we call "super-eccentric" Jupiters. Quantita- 
tive predictions and a strategy for detection using Kepler tar- 
gets are in §3. A brief summary is given in §4. 

2. STEADY-STATE DISTRIBUTION OF MIGRATING JUPITERS 

2.1. Basic assumptions 

Migration from large to small semi-major axis requires that 
energy is removed from the orbit. In HEM, tidal friction is re- 
sponsible for converting orbital energy into heat which is then 
radiated away from the system. In most cases, tidal friction 
in the planet removes energy much faster than tidal friction in 
the star. Since the planet's spin angular momentum is negli- 
gible compared to its orbital angular momentum, the orbital 
angular momentum per unit mass J is conserved during HEM 
(however, see discussion in §2.4), i.e.. 



J^ = G {M^,+Mp) fl (1 -e^) = est. 



(1) 



where a, Mp, M^, and G are the semi-major axis, planet mass, 
host star mass and gravitational constant, respectively. Thus 



[l-e ) = q{l +e) = flp = est. 



(2) 



where q is the periastron distance and flp is the final semi- 
major axis that the planet reaches when the eccentricity has 
decayed to zero. 

Let X = e or fl, and let Afj(X) dJ be the number of migrating 
planets with eccentricity or semi-major axis less than X and 
angular momentum in the interval {J, J+ dJ). We assume that 
all planets with eccentricity greater than some reference value 
Ci-ef are still migrating, and set A/}(ei-ef) = 0. We have argued 
that the distribution of migrating planets is in steady state and 
that orbital angular momentum is conserved during migration. 
Then the continuity equation requires 



X 



■dj\fj 



dX 



= Sj 



(3) 



where SjdJ is the current of migrating planets with angular 
momentum in the interval {J,J + dJ). This current is deter- 
mined by the properties of the source of highly eccentric long- 
period gas giants, which is assumed to be far {a^ I AU) from 
the region of phase space under consideration. 
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2.2. Orbital evolution: approximate treatment at high 
eccentricity 

We now describe an approximate analytic treatment of the 
orbital evolution and steady-state distribution at high eccen- 
tricity. For high eccentricity the shape of the orbit near peri- 
astron and the energy loss per periastron passage AE are both 
independent of e. Thus the orbit-averaged energy loss rate is 

. AE 1 

E= 

P 

where P is the orbital period. Since E ocl/a 
da .„ 



oc — oc a 
P 



-3/2 



(4) 



dt 



oc a 



1/2 



In the region of (e,a) space that contains a steady-state 
distribution of migrating planets on high-eccentricity orbits 
(q = a{l-e) < IQRq for Sun-like host stars) the number of 
migrating Jupiters per unit semi-major axis is found with the 
help of equation (3), 

dAfj est. _i/2 dAfj 



or 



,1/2 



da a dloga 
2.3. Orbital evolution: exact treatment 



(6) 



In order to study orbital evolution at small or moderate ec- 
centricity, some understanding of tidal dissipation is required. 
Unfortunately there is no robust theory of tidal dissipation in 
gas-giant planets, due both to the sparseness of observational 
calibration (only Jupiter and Saturn) and to theoretical diffi- 
culties in studying such weak dissipation (e.g., tidal Q ^ IQr' 
for the Jupiter-Io system). 

For illustration, we shall use the phenomenological ap- 
proach of Hut (1981), which follows Darwin in assuming that 
the tides lag their equilibrium value by a constant time r. By 
assuming pseudo-synchronous rotation (Hut's eq. 45) we find 
that the orbital evolution for a single planet is described by 



de 1 
— = — e(l ■ 
dt 2 



which is equivalent to 



— -- fl'/^ 
dt 



(7) 



(8) 



Here 
0.5 is 



where a = a/a^ = (1 - e^) ' and t = tjtu- 

is a dissipation time , k 
the planet's Love number, and Rp is its radius. This result 
assumes Mp <C M*. The function §(e) is given by 



^(e) = 



7+fg^ + 56g^+^g^ + ^g^ + ^g 
3(l+3e2+3g4) 

2.33 + 6. 12e^ 



(9) 



where the approximation in the final equation is accurate to 
better than 0.5% for all eccentricities between and 1. 

Equation (3) then implies that the number of planets per 
unit interval in angular momentum is given by 

dMAe) S] 



de 



\de/dt\ 



(10) 



For planets of a fixed density and a range of radii, the dissipation time 
scales as r (X l/{lcRp). Thus smaller planets have larger dissipation times. 
Nevertheless, we expect the dissipation to occur mostly in the planet rather 
than the star because the Love number Ic is much smaller in stars than planets. 




Fig. 1. — Cumulative distribution of migrating planets Afjie) with eccen- 
tricity e > eief = 0.1, from equation (11). The planets migrate along a track 
of constant orbital angular momentum /. The vertical axis also represents the 
time required for the eccentricity to decay to e^f . The magenta dot represents 
the current position of HD 80606b (e = 0.94, a = 0.45 AU) and the orange dot 
represents a hypothetical planet with semi-major axis a = 5 AU flowing along 
the same angular-momentum track as HD 80606b. The number of objects 
in the range e = 0.94-0.995 is comparable to the number of objects in the 
range e = 0.2-0.94. The plot is shown in normalized units so the curves ai'e 
independent of /. 
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Fig. 2. — The blue curve shows the distribution dj\fj / d\oga of mi- 
grating Jupiters per unit interval in angular momentum J from equation 
(10) and the dashed red curve shows the high-eccentricity approximation 
dj\f /dloga Qc a'/^ (eq. 6). The high-eccentricity approximation is accurate 
to ~ 20% for a/flp > 10 or e > 0.95. The magenta and orange points have 
the same meaning as in Figure 1 . The plot is shown in normalized units so 
the curves are independent of /. 



and 



Mj{e)= / dN-j(e) = Sj[t{e)-t(e,,f)] 



(11) 



where t(e) is the time required to migrate from some initial 
eccentricity near unity to e. The cumulative distribution Afj(e) 
for a single track in J, as determined from equations (7) and 
(11), is shown in Figure 1. 

Figure 2 displays the expected number of migrating Jupiters 
per unit log a, as obtained from equation (8), as well as the 
high-eccentricity approximation (6). For a fixed interval in 
orbital angular momentum, the number of migrating Jupiters 
is an increasing function of log a above e ~ 0.9. 
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2.4. The approximation of constant orbital angular 
momentum 

In HEM scenarios, gas giants are assumed to be born on 
nearly circular orbits and then acquire a large eccentricity af- 
ter exchanging their angular momentum with other planets or 
distant stellar companions, through close encounters in the 
former case or Kozai-Lidov (KL) oscillations in the latter. 
Therefore, orbital angular momentum is not a constant during 
the process of eccentricity excitation. Our analysis assumes 
that eccentricity excitation takes place at large semi-major 
axes (say, a > 5-10 AU) and focuses on the region of {a,e) 
space where substantial orbital decay has already occurred but 
the eccentricity is still moderate to large (say a < I AU and 
e > 0.2). It is not clear whether or not the approximation that 
migration takes place at constant angular momentum is accu- 
rate for all semi-major axes a < I AU. In what follows, we 
assess the validity of the constant J approximation of §§2.2 
and 2.3 in the presence of KL oscillations, which are the most 
likely cause of changes in the orbital angular momentum of 
the migrating planet. 

We performed many numerical integrations of the orbit- 
averaged restricted three body problem, including the effects 
of general relativity, tidal dissipation and tidal precession. 
Each simulation was initialized with a Jupiter-mass planet or- 
biting about a solar-mass star, placed in a nearly circular orbit 
with semi-major axis a ~ 3-5 AU. The system also contained 
a solar-mass companion star, placed at distances of 30- 1000 
AU with inclination of 85° < / < 90° relative to the planetary 
orbit. Only the quadrupole term of the companion's potential 
was considered. 

Typically, KL oscillations commenced at the start of the in- 
tegration, with large amplitude variations in orbital angular 
momentum /. Due to the strong dependence of tidal dissipa- 
tion on periastron distance q, dissipation takes place almost 
entirely in the vicinity of /niin, the minimum orbital angular 
momentum during a KL oscillation. As a result, the value of 
/min remains roughly fixed during migration. Precession due 
to general relativity acts to decrease the amplitude of the os- 
cillation in J (e.g., Blaes et al. 2002; Wu & Murray 2003; 
Fabry cky & Tremaine 2007). Once the oscillation amplitude 
in J is sufficiently small (< 10% in J), such that the dissipa- 
tion rate does not change considerably during each cycle, the 
mean value of / remains constant and equal to the final orbital 
angular momentum Jp. Therefore, the distribution of planets 
from then on can be computed by assuming a constant J = Jp. 
At this stage of migration, KL oscillations are considered to 
be "quenched." Quantitatively, KL oscillations are quenched 
at a semi-major axis aq given by 

- 1.8AU f y'V-V' 
^ V0.05AU/ \^ 0.4 J \MqJ 

[mqJ VlOOOAuJ [l-0.5^) ^^^^ 

where flpei and gpe, are the semi-major axis and eccentricity of 
the perturber and imin is the mutual inclination at the phase of 
the KL oscillation when J = J^m, while Mper is the perturber 
mass. For larger perturber mass or smaller semi-major axis, 
KL oscillations are quenched closer to the host star In partic- 
ular, nearby giant planets quench the oscillations at smaller 
radii than distant companion stars; for example, a Jupiter- 
mass perturber at 10 AU has ag « 0.2 AU. The major dif- 




FlG. 3. — Comparison of the analytic theory (red) of §2.3 with the migra- 
tion of HD 80606b, as depicted in Figure 1 of Fabrycky & Tremaine (2007). 
The blue lines show the density in eccentricity and log semi-major axis of 
a steady-state ensemble of planets that all follow the same trajectory as HD 
80606b. The planets only contribute to the density when their angular mo- 
mentum is small, in particular when a(l - e^) < 0.14,0. 11 and 0.08 AU. The 
normalization is chosen so that all the curves match at o/c/f = 5. Despite the 
presence of KL oscillations for a / > 30, the density is approximated by the 
analytic expressions derived in §2.3 to within a factor of two for all values of 
{e,a). 



ference between our various integrations of KL oscillations 
with tidal dissipation was the value of aq, due primarily to 
variations in the distance of the perturber 

During each integration we tracked the time that the planet 
spent in a bin of width Ay centered on the final angular mo- 
mentum Jp. We used this information to construct the eccen- 
tricity and semi-major axis distributions that would be present 
in this angular momentum bin in a steady-state population of 
planets following this migration path. We found that even 
in the presence of KL oscillations the constant J approxima- 
tion described in §§2.2 and 2.3 reproduced the distribution of 
planets to within a factor of two or better, so long as Ay was 
not more than about 20% of yp. That is, we found that the 
steady-state formulas derived in §2.3 by assuming constant y 
were still approximately valid, even though the orbital angu- 
lar momentum experiences large amplitude oscillations (see 
example of HD 80606b in §2.5). 

This surprising agreement results from the fact that tidal 
dissipation, and thus migration, occurs mostly when y ~ y^in, 
the minimum value of the angular momentum during a Kozai- 
Lidov cycle. As long as y^^in is close to yp, the final value 
of orbital angular momentum after quenching, then all of the 
migration takes place within the bin of width Ay. The time 
that the planets spend on Kozai-Lidov cycles outside the bin 
is irrelevant, since they do not migrate there. 

2.5. An example 

Consider the migration of the gas-giant planet HD 80606b, 
which currently has a = 0.45 AU and e = 0.93, correspond- 
ing to flp = a{\-e^) = 0.06 AU. The migration track has been 
modeled by Wu & Murray (2003) and Fabrycky & Tremaine 
(2007), who start with an initially nearly circular orbit with 
a = 5 AU, similar to Jupiter KL oscillations are excited by 
a distant companion star HD 80607. For the first Gyr of the 
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Fig. 4. — In solid blue are lines of constant angular momentum / that cor- 
respond to Pp = P(\ -e^fl^- = 1.25, 2.5, 5.0, 10.0, and20.0 days. Horizontal 
dotted lines are fixed values of 1/ 1 - e, separated by factors of two. The cor- 
responding values of e (0.2, 0.6. etc.,) are on the right hand side of the figure. 
The relative number of planets expected in each interval along a track in con- 
stant J is given by the numbers on the left (2.9, 1.0, etc.,). Tidal dissipation 
drives the planets to flow from the upper right comer to the lower left comer, 
along the lines of constant /. Black dots are cuiTent RV observations with 
Msin/ > 0.25 taken from exoplanet.org with the exception of HD 20782b 
(green dot) for which the updated e = 0.97 (O'Toole et al. 2009) is used. The 
magenta dot highlights the cuiTent position of HD 80606b and the orange dot 
represents a hypothetical planet at 5 AU flowing along the same track in J as 
HD 80606b. 



evolution the eccentricity oscillates between e^ax = 0.993 and 
Smin = 0.04-0.25. The amplitude of the KL oscillations then 
gradually decays; the oscillations are quenched by 2.8 Gyr, 
when the eccentricity is 0.97 and the semi-major axis is 2 AU 
- in agreement with equation (12) - and thereafter the eccen- 
tricity and semi-major axis decay at constant angular momen- 
tum, reaching zero eccentricity after 4 Gyr at a semi-major 
axis flF = 0.071 AU. 

In Figures 1 and 2 the magenta and orange points repre- 
sent the current position of HD 80606b and its hypothetical 
Jupiter-like "progenitor," respectively. 

Figure 3 shows in blue the density of planets in eccentric- 
ity and semi-major axis that would result from a steady-state 
ensemble of migrating planets with the same trajectory as HD 
80606b. Three plots are shown, for angular-momentum cut- 
offs y2 = G(M-,+Mp)a, with a, = 0.14,0. 11, and 0.08 AU (top 
to bottom). For comparison, the red lines show the analyti- 
cal predictions of §2.3. In the latter stages of migration, after 
the KL oscillations have been damped, the density matches 
the analytical estimate extremely well — this is not surprising 
since the assumption of evolution at constant angular momen- 
tum is satisfied to high accuracy. At larger eccentricities and 
semi-major axes, when KL oscillations are present, the blue 
curves are displaced from extrapolation of these theoretical 
predictions by up to a factor of two or so, but their shapes re- 
main similar as expected from the arguments of the preceding 
subsection. 

All three of Figures 1, 2 and 3 show that in a steady state, 
an unbiased sample of exoplanets containing one HD 80606b 
should contain more than one migrating planet with a sim- 
ilar periastron distance and even larger semi-major axis and 
eccentricity. 

3. OBSERVATIONS, PREDICTIONS, AND DISCUSSION 
3.1. Current observations 



We compile a list of all known exoplanets with Msin; > 
0.25My, of which the radial-velocity planets are displayed 
in Figure 4. The planetary parameters are taken from 
exoplanets . org with the exception of HD 20782b, the 
green dot, whose eccentricity was recently revised to e = 0.97 
(O'Toole et al. 2009). For each planet we compute = 
P{1- e^)^/^, the final orbital period that a planet would reach 
if its eccentricity decayed to zero at constant angular momen- 
tum. For reference a planet with Jupiter's orbital period and 

= lOd would have e = 0.991. 

The blue lines in Figure 4 are lines of constant orbital angu- 
lar momentum, along which planets flow from long to short 
orbital periods. The relative number of planets expected in 
each interval along a track in constant J is given by the num- 
bers on the left (2.9, 1.0, etc.,). If for example, one gas giant 
planet is found migrating in the Pp = 5-10 d bin within an ec- 
centricity range 0.9 < e < 0.95 (such as HD 80606b), then 
there should be > 1 planet migrating in this 5-10 d bin within 
an eccentricity range 0.95 < e < 0.975 as well. 

Among the known gas-giant exoplanets there is a signif- 
icant excess population having eccentricity consistent with 
zero and P = P^ < lOd, corresponding to a < 0.09 AU for a 
solar-mass host star In the HEM scenario, these are planets 
that were formed at several AU, excited to high eccentricity, 
migrated due to tidal friction, and have now completed the 
migration process. There is no such excess for larger periods; 
in HEM models this implies that tidal dissipation is unimpor- 
tant for planets with > 10 d and we discard these from our 
sample. 

From the remaining sample we calculate the number of 
"moderately eccentric" planets, which we define to be those 
with 0.2 < e < 0.6, and the number of "super-eccentric" plan- 
ets (e > 0.9). The first two lines of Table 1 summarize the 
number of gas-giant planets with moderate eccentricity, as 
detected by radial-velocity (RV) surveys and transit photom- 
etry with spectroscopic follow-up ("Transit-i-RV")^. Both the 
RV and Transit-i-RV categories yield a fraction of moderate- 
eccentricity planets that is roughly 1 /2 in the 5-10 d bin and 
much smaller, < 1/15, in the 3-5 d bin. The sharp decline 
for smaller values of P^ is consistent with the expectation that 
tidal dissipation is stronger for orbits with smaller periastron, 
so that in a steady state the fraction of planets in the migration 
pipeline is smaller. 

In HEM models both the moderately eccentric and super- 
eccentric planets are in steady-state migration and therefore 
the population ratio in these groups can be calculated using 
the models of §2.3. Thus we can predict the number of super- 
eccentric planets that should have been found in RV surveys; 
this prediction is shown in boldface in the third line along 
with the number of super-eccentric planets actually found so 
far in these surveys. Furthermore, the fractions of moderately 
eccentric RV and Transit-i-RV planets can be used with our 
models to predict the number of super-eccentric planets in the 

^ The planets in the "Transit+RV" line in Table 1 are obtained by using the 
following search string in exoplanets . org (Wright et al. 2011): 

MSINI [mjupiter] > 0.25 and PER[day] * (1 -ECC^)^/- > 3 

andPER[day]*(l-ECC-)'/2 < 10 and DISCMETH == "TRANSIT" 

They mostly consist of planets discovered in ground-based transit surveys 
(over 70%) as as well as a handful of objects discovered by the COROT and 
Kepler space-based telescopes and by RV surveys. All objects in this cate- 
gory have eccentricities that have been determined by spectroscopy. The vast 
majority of Kepler planets have no spectroscopic follow-up and hence are not 
included in this line. 
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TABLE 1 

Expected number of super-eccentric planets " 





3-5 days 


5-10 days 


RV (moderate/total) 


0/13 


4/9 


Transit+RV (moderate/total) 


3/46 


3/8 


RV (super-eccentric, theory/observed) 


vs. 


2-3 vs. 2 


Kepler (super-eccentric, theory) 


2 


3-5 



^"Moderate" denotes the eccentricity range 0.2 < e < 0.6. "super- 
eccentric" denotes e > 0.9, and "total" is < f < 0.6. Numbers in boldface 
are predictions obtained from the number of moderate-eccentricity planets 
in the RV and Transit-t-RV categories by assuming that all giant planets at 
small periods are formed by high-eccentricity migration (HEM) and apply- 
ing the model of §2.3. The period intervals (e.g., 3-5 days) refer to the 
final period - P{\ —e~)^^~, which is the period after HEM is complete 
and the orbit is circularized, assuming constant orbital angular momentum. 
Only planets with A/sin/ > 0.25 My are included in the statistics. The pre- 
dictions are for super-eccentric planets with orbital period P < 2 yr only. 
The results are based on queries to the exoplanets . org database in 
September 20 1 1 . The predictions do not account for eccentricity-dependent 
selection effects. 



Kepler sample (Borucki et al. 201 1); this prediction is shown 
in boldface in the last line of the Table. 

Before discussing these predictions we address selection 
effects. RV surveys may be biased against the detection of 
super-eccentric planets for at least two distinct reasons. First, 
sparse observations of high-eccentricity orbits are likely to 
miss the strong reflex velocity signal near periastron, lead- 
ing to non-detection of planets that would be detected at the 
same semi-major axis and smaller eccentricity, or to an un- 
derestimate of the eccentricity if the planet is detected (Gum- 
ming 2004, O'Toole et al. 2009). This bias only sets in for 
e > 0.6 so the fraction of moderately eccentric planets de- 
tected in RV surveys is much more reliable than the fraction 
of super-eccentric planets. Second, we are mostly concerned 
with avoiding biases against detecting highly eccentric plan- 
ets at a given angular momentum (i.e., along a given migration 
track), rather than at a given semi-major axis. Here there is 
an additional bias, since high-eccentricity orbits have longer 
periods and hence a periodic signal is harder to detect and 
characterize in a given time baseline. 

Despite these poorly understood selection effects, the pre- 
dicted and observed numbers of super-eccentric planets in RV 
surveys as shown in Table 1 are consistent. However, the 
numbers are too small to test the validity of HEM scenarios. 

For transit surveys the selection effects can be divided into 
geometric effects, which depend on the orientation of the ob- 
server relative to the star (i.e., whether or not a planet transits 
the star), and survey effects, which depend on properties of 
the survey (time baseline, photometric accuracy, etc.). With 
adequate baseline and signal-to-noise ratio, the most impor- 
tant selection effect is geometrical: the probability that the 
planet will transit is given by 

where and r are the stellar radius and heliocentric distance 
of the planet during transit and is an average over the az- 
imuth of the sightline, which is equal to an average over the 
true anomaly. Therefore, on a migration track of constant an- 
gular momentum, the geometric selection effects are indepen- 
dent of eccentricity. That is, a hot Jupiter progenitor with say, 
a = 1 AU and e = 0.975 has the same detection probability as a 
circularized hot Jupiter with a = 0.05 and e = 0. Survey selec- 



tion effects, in contrast, are biased against high-eccentricity 
orbits because the period is longer so there are fewer transits 
in a given period, and because the transits are shorter so the 
signal/noise ratio is smaller. However, these selection effects 
can be calculated and corrected for using the methods outlined 
in Borucki et al. (201 1), and should be relatively small since 
our sample is restricted to giant planets, which are relatively 
easy to detect. 

3.2. Predictions for Kepler 

Kepler observations yield the planetary radius R and orbital 
period P. To compare these results to our models we assume 
that our mass limit, 0.25My, corresponds to R = S/?® and that 
the planet population within a range of period is roughly the 
same as the population within the same range of P^ = P{1 - 
g2-)3/2 gjjj(;g most planets have small eccentricities. The most 
recent /Tep/er catalog (Borucki et al. 2011) contains 30 planets 
with R > 8/?9 and 3 d< P < 5d, and 16 in the same radius 
range with5 d<P < lOd.*" 

These can be combined with the results from ground-based 
surveys (line 2 of Table 1) to predict the number of moderate- 
eccentricity planets in each period range, and these numbers 
are combined with the steady-state HEM models in §2.3 to 
predict the numbers of super-eccentric planets in the Kepler 
catalog (line 4 of Table 1). These predictions should be under- 
estimates since Kepler will detect planets with longer periods 
as the mission progresses (the automated pipeline in Borucki 
et al. 201 1 only finds objects with P < 93 d). 

These results imply that Kepler should detect several super- 
eccentric (e > 0.9) giant planets (R > 8/?®) with orbital period 
<2 yr If an extended Kepler mission permits detections of 
planets with longer periods the predicted number is higher. A 
significant fraction of these could have e > 0.94 i.e., more ec- 
centric than HD 80606b, the current confirmed record-holder 

A typical member of this population, with a^ = Q.l AU and 
M > 0.25Mj, produces a stellar reflex velocity > 50m s"' near 
periastron. For objects on highly eccentric orbits with random 
orientations, most transits occur near periastron, where the re- 
flex velocity is close to the periastron value — quantitatively, 
over half of all transits occur when the reflex velocity is 
within 10% of the periastron velocity. Thus relatively few 
low-exposure RV measurements near the transit epoch should 
be sufficient to detect and measure a large eccentricity. We 
suggest that all of the Kepler gas-giant planetary candidates 
with periods above > 20 d (Borucki et al. 2011 list 34 ob- 
jects with R > 8/?® and periods between 20 d and 93 d) be 
followed spectroscopically near transit (with one or two addi- 
tional measurements at other phases to determine the systemic 
velocity). 

4. SUMMARY AND DISCUSSION 

The main result of this paper is that if hot Jupiters are 
formed by high-eccentricity migration (HEM), then there 
must be a steady-state current or flow of gas-giant planets mi- 
grating from large to small orbital periods. Since tidal dissipa- 
tion is required for HEM and is only effective out to distances 
of a few stellar radii in typical exoplanet systems, the cur- 
rent must consist of planets that either have periastrons of a 

* Note that the ratio of planets in these two period bins, 16/30 = 0.5, is 
larger than the corresponding ratio for ground-based surveys, 8/46 = 0.2 (the 
numbers are the same whether we use P or Pp). This result suggests that Ke- 
pler has less selection bias against long-period gas giants than ground-based 
surveys, which favors the detection of super-eccentric migrating planets. 
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few stellar radii or undergo Kozai-Lidov oscillations or other 
dynamical processes that regularly bring their periastrons to 
these small values. Moreover, because energy loss from tidal 
dissipation only occurs near periastron, the rate of energy 
loss on high-eccentricity orbits varies inversely with the or- 
bital period; thus for every migrating planet on a moderate- 
eccentricity orbit there should be many super-eccentric plan- 
ets (e > 0.9). We have computed the expected eccentricity 
and semi-major axis distribution of the steady-state current of 
migrating planets using Hut's (1981) model of tidal dissipa- 
tion and assuming pseudo-synchronous planetary spin. Our 
results indicate that several super-eccentric gas-giant planets 
should be present in the Kepler exoplanet catalog. These can 
be discovered, if present, by a program of radial-velocity mea- 
surements on the Kepler planets with the largest diameters and 
the longest periods. 

The absence of a significant number of super-eccentric mi- 
grating Jupiters in this sample would imply either that HEM 
is not an ingredient of the formation process for most hot 
Jupiters, or that our migration model is oversimplified. In 
particular, we assume that migration occurs at constant or- 
bital angular momentum but argue that our results should be 
approximately correct even in the presence of Kozai-Lidov 



oscillations or other processes. 

The simple HEM model described here, whose central com- 
ponents are the steady-state approximation and the assump- 
tion that migration occurs at constant angular momentum, 
provides a preliminary framework for the exploration of the 
dynamics of HEM. A thorough exploration of this dynamics 
should establish whether our simplified model is accurate and 
enable a definitive observational test of whether hot Jupiters 
form through HEM. 
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